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^ I 

■ Abstract: This paper is devoted to study a model of individual clustering with two 

, specific reproduction rates in two space dimensions. Given q > 2 and an initial condition 

I in W^''^{Q), the local existence and uniqueness of solution have been shown in [6]. In this 

' paper we give a detailed proof of existence of global solution. 

, , 1 Introduction 
CLh 

I In the present work, we deal with a model of individual dispersing of individual with an 

r~| ' additional aggregation mechanism introduced in [5]. Given a sufficiently smooth function 

. E, parameters (5 G (0, 1), e > and r > 0, the equations take the form 

dtu = 6 Au-V ■ {u u) +r u E{u), xeQ,t>0 , , 

-eAu + u = VE{u), xen,t>0 

Q ! in an open bounded domain C M^, where u{t,x) > 0, u{t,x) £ E? and E denote the 

G> population density, the average velocity of dispersing individuals, and the individual net 

reproduction rate, respectively. In this model, the individuals are assumed to disperse ran- 
domly in space {6 An) with a bias —V • {u u) in the direction of increasing reproduction 



o 



\ rate, the term e Au acting as a moUifier to smooth out any sharp local variation in \/E{u). 



In [5], we supplement (_1) with no-flux boundary conditions 

dnU = n ■ u = 0, X e (90, t > 0, 



X 

. where n is the outward unit normal of dQ. and dnU = n • Vu. However, to guarantee the 

well-posedness of the elliptic system for u: in two space dimensions, we should append the 
following condition given in [3, 4, 7] 

drrU X n = 0, X £dn, t>0, (3) 

where dnu: = {dnUJi, 9n'^2) = • VtJi, n • 'VUJ2) for the vector field u = (wi, W2)- in other 
words, (3) means that dnUJ is parallel to n, where v x u = vi U2 — ui V2- 

Given q > 2, and an initial condition uq G W^''^{Q,), the existence and uniqueness of 
a nonnegative and maximal solution of {!), (2) and (3) have been shown in [6], and the 
purpose of this paper is to prove the global existence of solution when E{u) has the two 
specific forms suggested in [5], namely 

E{u) = (1 - u) {u- a) (4) 
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for some a G (0, 1), or 



E{u) = l-u. 



(5) 



For these choices of reproduction rates, global existence has been shown in [6] in one space 
dimension and the purpose of this work is to prove that the solutions are global as well in 
two space dimensions. As in the one-dimensional case, the starting point of the analysis 
is an L°°(L^) estimate on u and an estimate on V • a;. Combining the latter with 
Gagliardo-Nirenberg inequality gives L°^{LP) estimates on u for any p > 2. This then 
allow us to obtain an L°°{Lp') bound on Vu which in turn gives an L°° bound on u by 
elliptic regularity. 

The paper is organized as follows. In section 2, we state the global existence results, 
and focus on the two specific forms of E: The "bistable case" (4) see Theorem 2.2, and 
the "monostable case" (5), see Theorem 2.3. In section 3, we recall the local existence 
result obtained in [6] and we give some properties of the elliptic system for In section 
4, we turn to the global existence issue in the bistable case. The proof starts from the 
L°°(L^) estimate for u, an estimate for u and an L? estimate on V • cj obtained from 
a suitable cancellation between the coupling terms in the u and equations, then, for 
p > 2, we derive an L°°{LP) estimate for u. Then we use Lemma A.l of [8] to derive an 
L°° estimate of u and we end the proof by an L°°{L'^) estimate on Vu. This ensures global 
existence. In section 5, we prove the global existence in the monostable case. The proof 
is quite similar to that of the previous case, except for the first estimate. 



2 Main result 

We first define the notion of solution to (l)-(3) to be used in this paper. 

Definition 2.1. Let T > 0, q > 2, and an initial condition uq G W^''^{Q) . A strong 
solution of (JL}-(3) on [0,T) is a function 

u€C{[o, T), w^''^{n)) n c ((o, r), w^'''{Q)) , 



such that 




6 Au-V ■ {u u) + r u E{u), a.e. in [0, T) x 

lio(x), a.e. in Q. (6) 

0, a.e. on [0,r) x dQ, 



where, for all t G [0, T), u:{t) is the unique solution in W'^''^{Q.) of 

-eAu{t)+u{t) = VE{u{t)) a.e. in 17 

ijj{t) • n = dnU>{t) X n = a.e. on dO. 



(7) 



In the following theorem we give the global existence of solution to (l)-(3) in the 
bistable case, that is when E{u) = (1 — u){u — a), for some a G (0, 1). 

Theorem 2.2. Let q > 2, and assume that uq is a nonnegative function in W'^''^{Q), 
and E{u) = {1 — u){u — a) for some a G (0, 1). Then (JL)-(3.) has a global nonnegative 
solution u in the sense of Definition 2.1 . 

The proof starts with a suitable cancellation of the coupling terms in the two equations 
which gives an estimate for u in L°°{L'^) and for V-o; in L^. Using the Gagliardo-Nirenberg 
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inequality (13) we derive, for p > 2 an L°°{LP) estimate for u. Then by the regularity 
properties of the second equation in (1) we obtain an L°° bound of lj. Combining these 
estimates and Lemma A.l of [8] provide us with an L°° estimate for u which is used to 
show an L°°{L'^) estimate for Vn. This proves that the solution cannot explode in finite 
time. 



Next, we turn to the global existence issue in the monostable case, that is when E{u) = 
l-u. 

Theorem 2.3. Let q > 2, and assume that uq is a nonnegative function in W^''^{Q), 
and E{u) = [\ — u). Then (l)-(3) has a global nonnegative solution u in the sense of 
Definition 2.1. 

The proof of the previous theorem follows the same lines as that of Theorem 2.2. As 
in the bistable case, there is a cancellation between the two equations which provide us an 
L°°{L\ogL) bound on u and an bound for V • tJ as a starting point. 



3 Well-posedness 

Throughout this paper and unless otherwise stated, we assume that 

5 £ (0,1), e>0, r > 0. 
We first recall some properties of the strong solution of the following system, 

—e A(jL> + (jL> = f, in Q, 

u n = 0, on dn, (8) 

dnu: X n = 0, on dO,, 

where / G (L^(r2))^ and p > 1. The strong solutions of (8) is solving (8) a.e. in Q,. In this 
direction the existence and uniqueness of the strong solution to (8) are proved in [7] : 

Theorem 3.1. For f € (LP(il))^ with 1 < p < oo, (8) has a unique solution in (H^^'^'($7))^ 
such that 

\Mw-,. < ^ WfWp, (9) 

where K{p) = K{p, Q). 

In other words, the strong solution has the same regularity as elliptic equations with 
classical boundary conditions. 

Thanks to [6] we recall the existence and uniqueness result of the maximal solution of 
(l)-(3). 

Theorem 3.2. We assume that E E (7^(M), let p > 2 and a nonnegative function 

uq G W^'^{^). Then, for some Tmax G (0,oo], there is a unique nonnegative maximal 

solution 

uec {[o,T^^^),w''P{n)) n c {{o,T„,,^),w^'P{n)) (lo) 

to (l)-(3) in the sense of Definition 2.1 . Moreover, if for each T > 0, there is CiT) such 
that 

\\u{t)\\wuv < C{T), for all t e [o,r] n [o,r,nax), 
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then Tmax = oo. In addition, u satisfies 

u{t, x) = (^e'^^ ^) uo(x)) + e(*-^)(^ ^) [- V ■ {u lj) + r u E{u)] (s, x) ds, (11) 

for {t,x) G [OjTjnax] X where (e* ^•') denotes the semigroup generated in LP{yL) by 
5 A with homogeneous Neumann boundary conditions. 

We recall that there is C > such that 

||e* ^'^v\\wi,v < C \\v\\wi-P, and ||Ve* ^h\[p < C 5'^ t'^ \\v\\p. (12) 
Also in several places we shall need the following Gagliardo-Nirenberg inequality 

\\u\\p < C ||n||^i,2 \\u\\\~\ with 9 = u G W^''^{n) (13) 

which holds for all p > 1 and q G [l,p]- Also we use the following singular Gronwall lemma 
(see [1, Theorem 3.3.1]). 

Lemma 3.3. Given a,/3 £ [0,1) , there exists a positive constant c := c{a,(3) such that 
the following is true: 

Iff : (0, T) — >M. satisfies 

't^t^ fit)] GLi-((0,T),M), (14) 

and ^ 

f{t)<At-^ + B [ .^—f{s)ds, a.a.te{0,T), (15) 

where A and B are positive constants, then f{t) < C{T), for all t G (0,T), where C 
depends only on T, a, f3, and 7. 

4 Global existence 

4.1 The bistable case: E{u) = (1 — u){u — a) 
We recall the system 



(16) 



dtu = 5 Au — V ■ [u uj) + r u {I — u) {u — a), x e 0,,t > 

-eAu + u = [-2 n + (a + 1)] Vn, xen,t>0 

dnU = , u> ■ n = dnU X n = 0, a; G dQ, t > 

u(0, x) = no(x), X G ri. 

for a some a G (0, 1), and uq G W^''^{Q) for some q> 2. 

Dmce E G C2(IR), Theorem 3.2 ensures that there is a maximal solution of (16) in 
C([0,T„,a.),W^''?(f^)) nC((0,T^ax),W'''(f^)) for q > 2. 

We begin the proof by the following lemmas which gives some estimates on u and u. 
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Lemma 4.1. Let the same assumptions as that of Theorem 2.2 hold, and u he the non- 
negative maximal solution of (16). Then for all T > there exists Ci{T) > 0, such that 
u and u satisfy the following estimates 



/ \\Vu{s)\\ids<Ci{T), for alHe [0,r]n[0,r^ax), (17) 
Jo 

and ^ 

[ {\\V-u{s)\\l + \\u{s)\\l) ds<Ci{T) foralHG[0,r]n[0,T^ax). (18) 
Jo 

Proof. We multiply the first equation in (16) hy 2 u and integrate it over 0, to obtain 

^ ( \u\'^ dx = -2 5 f \Vu\'^ dx + 2 [ uu-Vudx + 2r [ u"^ E{u) dx. (19) 
Jn Jq Jn Jn 

We multiply now the second equation in (16) by u and integrate it over Q. We note that 
the boundary conditions for u guarantee that u is tangent to 9$7 while dnCiJ is normal to 
dQ. Consequently, dnUJ • cj = on d^l and it follows from an integration by parts that 

—e / Au ■ u> dx = e / |V • cjp dx — e [ (Vwi ■ n) loi + {Vlo2 ■ n) L02] da 
Jn Jn Jan 



e I \V ■ dx. 
Jn 



We thus obtain 



e / |V • cj|^ + / |cj|^ dx = -2 u u ■ Vu dx + {a + 1) u ■ Vu dx. (20) 
Jn Jn Jn Jn 

At this point we notice that the cubic terms on the right hand side of (19) and (20) cancel 
one with the other, and summing (20) and (19) we obtain 

+ e ||V • cj||2 + 11^112 + 2 5 \\Vu\\l = 2 r u^ E{u) dx + {a + 1) / cj-Vu dx. 

Jn Jn 



l|2 

di" 



We integrate by parts and use Cauchy-Schwarz inequality to obtain 

(a + 1) [ u - Vu dx = -(a + 1) / u V ■ u dx < ll-uHo + - ||V • uj\\l. 

^ Jn Jn - 2 e " 2 

On the other hand, E{u) < if u ^ (a, 1) so that 

/ u"^ E{u) dx < \Vt\ (1 - a). 
Jn 

The previous inequalities give 

-^llullo + - IIV • culln + llculli + 2 5 llVnllo < llnlln +2 101 r (1 - a), 

dt 2 2 e 

Therefore, for all T > there exists Ci(T) such that (17) and (18) hold. □ 

Lemma 4.2. Let the same assumptions as that of Theorem 2.2 hold, and u he the non- 
negative maximal solution of (16). Then for allT > there exists C2{T,p) > 0, such that 
for p >2 

\\u{t)\\p<C2{T,p) for aUtG [o,r]n[o,r^ax), (21) 

['\\VuHs)\\l ds < C2{T,p) for all t G [0,r] n [0,r,„ax). (22) 
Jo 
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Proof. We multiply the first equation in (16) hy p ^, integrate with respect to x, and 
integrate by parts. The boundary terms vanish and we obtain 

dt" - p ' Jn 

+ r p u^^^ E{u) dx. 
Jn 

By Cauchy-Schwarz inequality we obtain 



^IHI^ < "'^|f"'^ ||Vni||i + (p-l)||V-c.|b||ni||l (23) 
+ r p {1 — a) 

Using the Gagliardo-Nirenberg inequality (13) we have 

||n2||4 < C ||'U2|||j,i,2 \\u2\q. (24) 
Substituting (24) in (23), and by Young inequality we obtain 

d — 4 6 (v — 1) p n p p 

-nlMl < ^ l|Vn2||^ + C (p-l) ||V-cj||2 ||m2||i4^i,2 ||m2||2 

dr ^ p 

+ r p (l-a) \n\ 

p p 
+ C{p) \\um + C{p) 

< l|Vn§||i + C(p) \\u\\P + C{p) ||V-^||i ||«||^ + C(p). 

Next, integrating the above inequality in time, and using (18) yield that there exists 
C2{T,p) such that (21) and (22) hold. □ 

Lemma 4.3. Let the same assumptions as that of Theorem 2.2 hold, and u be the non- 
negative maximal solution of (16). Then for all T > there exists C^{T) > 0, such 
that 

||Vn(t)||2 < C3(T) for all t G [0,r] n [0,T^ax). (25) 

Proof. We multiply the first equation in (16) by —Au, integrate over Q, and use Cauchy- 
Schwarz, and Young inequalities and (21) to obtain 



- — \ 
2dt' 



Vulll = -^||A„||i+/(V„.. + V..„) AW.-./„(l-„)(„-„)AW. 

Jn Jn 

\\Au\\l + C [ |Vu|2 dx 
Jn 



< -5 \\Au\\l 



+ C \V ■ uj\^ \u\^ dx + C r \\u {I - u) {u - a)\\l + -\\Au\\l 
Jn 4 

[ |VnHa;|2dx + C / \V ■ u:\'^ \u\'^ dx + C (T) . (26) 
Jn Jn 



6_ 
4 



< -^\\Au\\l + C 



To go further requires to improve the estimate on u and V • cj. For that purpose, we use 
Lemma 4.1 and Lemma 4.2 for p = 4 to obtain for all T > 
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/ \\VE{u)is)\\lds<{a + lf [ \\Vu{s)\\l ds + 4 [ \\u Vu{s)\\l ds < C2{T), (27) 
Jo Jo Jo 

for all t G [0,T] (1 [0,Tmax)- Consequently, VE{u) is bounded in L^iiO,t) x O). By 
Theorem 3.1 and the continuous embedding of VF^'^(r2) in W^''^{Q), and W^''^{Q) in L°°{Q), 
we have 

Halloo + ||V • CJ||4 < C IMwlA < C \\u{s)\\w2,2 < C ||V^(U)||2, 

which together with (27) implies that 

[\\\u^{s)\\l + \\V-u{s)\\l) ds<C f\\VE{u){s)\\lds<C2{T), (28) 
Jo JO 

for all t G [0,T] n [0,Tmax)- Using Holder Inequality , (26) becomes 

C II^IlL / iV-up dx + C ||V • ^lll \\u\\l + C{T). 
Jn 



1 d ,,9 6 ,, , ,,9 

--r 9 < -T Ati 9 
2dt" - 4 " 



Next we integrate the above inequality in time, and use (21) for p = 4, and (28) to obtain 
\\Vu{t)\\l < WVuoWl + C t\\u{s)\\l \\Vn{s)\\l ds + C{T), 



using (28) again, we have thus proved (25). □ 

Lemma 4.4. Let the same assumptions as that of Theorem 2.2 hold, and u he the non- 
negative maximal solution of (16). Then for all T > there exists C4(T) > 0, such 
that 

\\^{t)\\oo < C4{T), for ah t G [o,r] n [o,r^ax). (29) 

Proof. It follows from (21), (25) and Holder inequality, that there exists C(T) > such 
that 

\\(-2u + a + l) Vnlh < C(T) \\ - 2 u + a + llL llVulb < C(T). 

2 

Consequently, Theorem 3^1. ensures that 

ll^ll 2 3 < C ||[-2 u + a + 1] Vu\\3 < C{T), for ah t G [0,T] n [0,rmax)- 

Using the continuous embedding of (0) in L°^{Q) we have thus proved (_29). □ 

Next, thanks to lemma A.l in [8] we can derive a uniform bound for u. 

Lemma 4.5. Let the same assumptions as that of Theorem 2.2 hold, and u be the non- 
negative maximal solution of (16). Then for all T > there exists C^{T) > 0, such 
that 

\\u{t)\\oo < C5{T), for ah t G [o,r] n [o,r^ax). (30) 
Proof. We can see that the function u solves 

dtu = V-{5Vu) + V-f + g, a.e. in [0,Ti„ax) X J] 
u{0,x) = uo{x), a.e. in (31) 

dnU = 0, a.e. on [0, Tmax) x dil, 



7 



where / = 



—u ijj and g = r u E{u). 



By the regularity (10) of u, and by the continuous embeddings of W^'P{Q) in C{Q) 
and of W'^'P{n) in C{n) for p > 2 we obtain that f = -u u e C {{0,T^a^); C(0)) and 
V/ G C ((0,rmax); C(0)). Using (10) and the continuous embeddings of W^^p{Q) in C{Q) 
for p > 2 again, we see that g = r u E{u) G C ((0,Tinax) x On the other hand, 
f ■ n = —u u ■ n = on dQ x (0, Tmax)- 

Thanks to (21), we have 

Mmpo < c{T), for all t G [o,r] n [o,r^ax) 

where po = 9, while (21) and (29) yield 

\\f{t)\\g, < C{T) and \\g\\g, < C{T), for ah t G [0,T] n [0,r^ax) 

where (?i > 4 and q2 = 3 > 2. Since po > 1 — ^q^-s = 9 > 0, we can apply lemma 

A.l in [8], and the estimate (30) holds. □ 

Lemma 4.6. Let the same assumptions as that of Theorem 2.2 hold, and u be the non- 
negative maximal solution of (16). Then for all T > there exists Cq{T) > 0, such 
that 

\\yu{t)\\g < Ce{T), for all t G [0,T] n [0,r^ax). (32) 
Proof Using (11), (12), (30) and (29) we have for g > 2 

\\Vu{t)\\g < C\\uo\\wi,,+Ci [ it - s)~^ \\V ■ {u u){s)\\q ds 

Jo 

+ rCi [ \\V{u E{u))is)\\g ds 







r 1 

< C \\uo\\wi,q + Ci (t - sy^ \\u{s)\\oo \\V ■ u;{s)\\g ds 

Jo 

+ Ci {t-s)-^\\u{s)\\oo\\Vuis)\\qds (33) 
Jo 

+ rCi [ \\{-3u^ + 2{a + l)u-a)is)\\oo\\Vu{s)\\gds. 
Jo 

< C{T) (^1 + jjt- s)-5(||V • u{s)\\g + ||Vn(s)||,) ds + ||Vn(s)||, ds^ . 

By Theorem 3.1 we have 

\\V-Lj{s)\\g < \\u:{s)\\w2,, < K{q) ||(-2n + a + l)(s)||oo \\yu{s)\\g. (34) 
Substituting (34) into (33) and using (30) we obtain 

||Vn(i)||g < ^(r) {l + j\t - s)-5 \\Vu{s)\\g ds + \\Vu{s)\\g ds^ . 

Now, we obtain (32) by applying Lemma 3.3 with /3 = 0, and a = ^. □ 
Thanks to these lemmas we can now prove the main Theorem 2.2 
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Proof of Theorem 2.2. For all T > 0, Lemma 4.2 and Lemma 4.5 ensure that for g > 2 

< C(r), for all t E [0, T] n \^,T,,,^^) 

which guarantees that u cannot explode in VF^''?(r2) in finite time and thus 

that Tmax = oo. □ 

4.2 Monostable case 

For this choice of i?, system (l)-(3) now reads 



' dtu = (5 Au - V • (u cj) +r u (1 - n), 3;eO,t>0 

-eA^ + cj = - Vn, xGO,t>0 

(9n^ = , a; • n = Srjo; x n = 0, x G 5r2, t > 

n(0, x) = iio(a;), X G 0, 



(35) 



when G W^^'^iVC) for some q > 2. Since E G C^(M), Theorem 3.2 ensures that there is 
a maximal solution of (35) in C ([0,T„iax), W^^i{9.)) n C ((0, T^ax), ^^^'^(J^)) for g > 2. 

To prove Theorem 2.3 we need to prove the following lemma 

Lemma 4.7. Let the same assumptions as that of Theorem 2.3 hold, and let u he the 
maximal solution of (35) • Then for all T > 0, there exists Cj{T) > such that 

[\\\u{s)\\l + \\V •u{s)\\l) ds<C7iT), for ah t G [0,r]n[0,r^ax). (36) 
Jo 

Proof. The proof goes as follows. On the one hand, we multiply the first equation in (35) 
by (log n+ 1) and integrate it over Q, the boundary terms vanish. Since u {1 — u) log n < 
and u {1 — u) < 1, 

-r f u logu dx = — [ Vti — u u) ■ (— Vu) dx + r [ u (1 — u) (logu + 1) dx 
citjQ Jn u Jn 

< - f - |Vn|^ dx+ f u-Vu dx + \n\ r. (37) 



On the other hand, we multiply the second equation in (35) by u, integrate it over fi. We 
note that the boundary conditions for u guarantee that u is tangent to d^l while dnUJ is 
normal to dil,. Consequently, dnUJ • a; = on dfi and it follows from an integration by 
parts that 

e 1 \V dx+ I \u>\'^ dx = - u-Vu dx. (38) 
Jn Jn Jn 

Adding (38) and (3_7) yields 

^ [ u logu dx + e ||V -culll + <-A5 f |VVup dx + \n\ r. (39) 

dt Jn Jn 

Finally, (36) is obtained by a time integration of (39). □ 

Proof of Theorem 2.3 . Thanks to (36) , we now argue as in the proof of Lemma 4.2, Lemma 
4.3 . Lemma 4.4, Lemma 4.5 and Lemma 4.6 to get that for q > 2 

||n(t)||oo + ||Vu(t)||g < C(r), for ah t G [o,r] n [o,r^ax). 

Thus, the maximal solution u of (35) cannot explode in finite time. □ 
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